Transcriptomic Profiling of Reproductive Age Marmoset Monkey Ovaries.

The decline in ovarian reserve and the aging of the ovaries is a significant concern for women, particularly in the context of delayed reproduction. However, there are ethical limitations and challenges associated with conducting long-term studies to understand and manipulate the mechanisms that regulate ovarian aging in human. The marmoset monkey offers several advantages as a reproductive model, including a shorter gestation period and similar reproductive physiology to that of human. Additionally, they have a relatively long lifespan compared to other mammals, making them suitable for long-term studies. In this study, we focused on analyzing the structural characteristics of the marmoset ovary and studying the mRNA expression of 244 genes associated with ovarian aging. We obtained ovaries from marmosets at three different reproductive stages: pre-pubertal (1.5 months), reproductive (82 months), and menopausal (106 months) ovaries. The structural analyses revealed the presence of numerous mitochondria and lipid droplets in the marmoset ovaries. Many of the genes expressed in the ovaries were involved in multicellular organism development and transcriptional regulation. Additionally, we identified the expression of protein-binding genes. Within the expressed genes, VEGFA and MMP9 were found to be critical for regulating ovarian reserve. An intriguing finding of the study was the strong correlation between genes associated with female infertility and genes related to fibrosis and wound healing. The authors suggest that this correlation might be a result of the repeated rupture and subsequent healing processes occurring in the ovary due to the menstrual cycle, potentially leading to the indirect onset of fibrosis. The expression profile of ovarian aging-related gene set in the marmoset monkey ovaries highlight the need for further studies to explore the relationship between fibrosis, wound healing, and ovarian aging.

Changes in oral health during aging in a novel non-human primate model.

Oral health plays a significant role in the quality of life and overall well-being of the aging population. However, age-related changes in oral health are not well understood due to challenges with current animal models. In this study, we analyzed the oral health and microbiota of a short-lived non-human primate (i.e., marmoset), as a step towards establishing a surrogate for studying the changes that occur in oral health during human aging. We investigated the oral health of marmosets using cadaveric tissues in three different cohorts: young (aged ≤6 years), middle-aged, and older (>10 years) and assessed the gingival bacterial community using analyses of the V3-V4 variable region of 16S rRNA gene. The oldest cohort had a significantly higher number of dental caries, increased dental attrition/erosion, and deeper periodontal pocket depth scores. Oral microbiome analyses showed that older marmosets had a significantly greater abundance of Escherichia-Shigella and Propionibacterium, and a lower abundance of Agrobacterium/Rhizobium at the genus level. Alpha diversity of the microbiome between the three groups showed no significant differences; however, principal coordinate analysis and non-metric multidimensional scaling analysis revealed that samples from middle-aged and older marmosets were more closely clustered than the youngest cohort. In addition, linear discriminant analysis effect size (LEFSe) identified a higher abundance of Esherichia-Shigella as a potential pathogenic biomarker in older animals. Our findings confirm that changes in the oral microbiome are associated with a decline in oral health in aging marmosets. The current study suggests that the marmoset model recapitulates some of the changes in oral health associated with human aging and may provide opportunities for developing new preventive strategies or interventions which target these disease conditions.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Whole Genome Analysis of SNV and Indel Polymorphism in Common Marmosets (Callithrix jacchus).

The common marmoset (Callithrix jacchus) is one of the most widely used nonhuman primate models of human disease. Owing to limitations in sequencing technology, early genome assemblies of this species using short-read sequencing suffered from gaps. In addition, the genetic diversity of the species has not yet been adequately explored. Using long-read genome sequencing and expert annotation, we generated a high-quality genome resource creating a 2.898 Gb marmoset genome in which most of the euchromatin portion is assembled contiguously (contig N50 = 25.23 Mbp, scaffold N50 = 98.2 Mbp). We then performed whole genome sequencing on 84 marmosets sampling the genetic diversity from several marmoset research centers. We identified a total of 19.1 million single nucleotide variants (SNVs), of which 11.9 million can be reliably mapped to orthologous locations in the human genome. We also observed 2.8 million small insertion/deletion variants. This dataset includes an average of 5.4 million SNVs per marmoset individual and a total of 74,088 missense variants in protein-coding genes. Of the 4956 variants orthologous to human ClinVar SNVs (present in the same annotated gene and with the same functional consequence in marmoset and human), 27 have a clinical significance of pathogenic and/or likely pathogenic. This important marmoset genomic resource will help guide genetic analyses of natural variation, the discovery of spontaneous functional variation relevant to human disease models, and the development of genetically engineered marmoset disease models.

Extensive Independent Amplification of Platy-1 Retroposons in Tamarins, Genus Saguinus.

Platy-1 retroposons are short interspersed elements (SINEs) unique to platyrrhine primates. Discovered in the common marmoset (Callithrix jacchus) genome in 2016, these 100 bp mobile element insertions (MEIs) appeared to be novel drivers of platyrrhine evolution, with over 2200 full-length members across 62 different subfamilies, and strong evidence of ongoing proliferation in C. jacchus. Subsequent characterization of Platy-1 elements in Aotus, Saimiri and Cebus genera, suggested that the widespread mobilization detected in marmoset (family Callithrichidae) was perhaps an anomaly. Two additional Callithrichidae genomes are now available, a scaffold level genome assembly for Saguinus imperator (tamarin; SagImp_v1) and a chromosome-level assembly for Saguinus midas (Midas tamarin; ASM2_v1). Here, we report that each tamarin genome contains over 11,000 full-length Platy-1 insertions, about 1150 are shared by both Saguinus tamarins, 7511 are unique to S. imperator, and another 8187 are unique to S. midas. Roughly 325 are shared among the three callithrichids. We identified six new Platy-1 subfamilies derived from Platy-1-8, with the youngest new subfamily, Platy-1-8c_Saguinus, being the primary source of the Saguinus amplification burst. This constitutes the largest expansion of Platy-1 MEIs reported to date and the most extensive independent SINE amplification between two closely related species.

Population genetics of marmosets in Asian primate research centers and loci associated with epileptic risk revealed by whole-genome sequencing.

The common marmoset ( Callithrix jacchus) has emerged as a valuable nonhuman primate model in biomedical research with the recent release of high-quality reference genome assemblies. Epileptic marmosets have been independently reported in two Asian primate research centers. Nevertheless, the population genetics within these primate centers and the specific genetic variants associated with epilepsy in marmosets have not yet been elucidated. Here, we characterized the genetic relationships and risk variants for epilepsy in 41 samples from two epileptic marmoset pedigrees using whole-genome sequencing. We identified 14 558 184 single nucleotide polymorphisms (SNPs) from the 41 samples and found higher chimerism levels in blood samples than in fingernail samples. Genetic analysis showed fourth-degree of relatedness among marmosets at the primate centers. In addition, SNP and copy number variation (CNV) analyses suggested that the WW domain-containing oxidoreductase ( WWOX) and Tyrosine-protein phosphatase nonreceptor type 21 ( PTPN21) genes may be associated with epilepsy in marmosets. Notably, KCTD18-like gene deletion was more common in epileptic marmosets than control marmosets. This study provides valuable population genomic resources for marmosets in two Asian primate centers. Genetic analyses identified a reasonable breeding strategy for genetic diversity maintenance in the two centers, while the case-control study revealed potential risk genes/variants associated with epilepsy in marmosets.

Adeno-associated viral vectors for functional intravenous gene transfer throughout the non-human primate brain.

Crossing the blood-brain barrier in primates is a major obstacle for gene delivery to the brain. Adeno-associated viruses (AAVs) promise robust, non-invasive gene delivery from the bloodstream to the brain. However, unlike in rodents, few neurotropic AAVs efficiently cross the blood-brain barrier in non-human primates. Here we report on AAV.CAP-Mac, an engineered variant identified by screening in adult marmosets and newborn macaques, which has improved delivery efficiency in the brains of multiple non-human primate species: marmoset, rhesus macaque and green monkey. CAP-Mac is neuron biased in infant Old World primates, exhibits broad tropism in adult rhesus macaques and is vasculature biased in adult marmosets. We demonstrate applications of a single, intravenous dose of CAP-Mac to deliver functional GCaMP for ex vivo calcium imaging across multiple brain areas, or a cocktail of fluorescent reporters for Brainbow-like labelling throughout the macaque brain, circumventing the need for germline manipulations in Old World primates. As such, CAP-Mac is shown to have potential for non-invasive systemic gene transfer in the brains of non-human primates.

Optimization of piggyBac transposon-mediated gene transfer method in common marmoset embryos.

Generating non-human primate models of human diseases is important for the development of therapeutic strategies especially for neurodegenerative diseases. The common marmoset has attracted attention as a new experimental animal model, and many transgenic marmosets have been produced using lentiviral vector-mediated transgenesis. However, lentiviral vectors have a size limitation of up to 8 kb in length for transgene applications. Therefore, the present study aimed to optimize a piggyBac transposon-mediated gene transfer method in which transgenes longer than 8 kb were injected into the perivitelline space of marmoset embryos, followed by electroporation. We constructed a long piggyBac vector carrying the gene responsible for Alzheimer's disease. The optimal weight ratio of the piggyBac transgene vector to the piggyBac transposase mRNA was examined using mouse embryos. Transgene integration into the genome was confirmed in 70.7% of embryonic stem cells established from embryos injected with 1000 ng of transgene and transposase mRNA. Under these conditions, long transgenes were introduced into marmoset embryos. All embryos survived after transgene introduction treatment, and transgenes were detected in 70% of marmoset embryos. The transposon-mediated gene transfer method developed in this study can be applied to the genetic modification of non-human primates, as well as large animals.

The Verriest Lecture: Pathways to color in the eye and brain.

In common with the majority of New World monkeys, marmosets show polymorphic color vision by allelic variation of X-chromosome genes encoding opsin pigments in the medium/long wavelength range. Male marmosets are thus obligate dichromats ("red-green color blind"), whereas females carrying distinct alleles on X chromosomes show one of three trichromatic phenotypes. Marmosets thus represent a "natural knock-out" system enabling comparison of red-green color vision in dichromatic and trichromatic visual systems. Further, study of short-wave (blue) cone pathways in marmosets has provided insights into primitive visual pathways for depth perception and attention. These investigations represent a parallel line to clinical research on color vision defects that was pioneered in studies by Guy Verreist, whom we honor in this eponymous lecture.

Metagenomic analysis reveals novel dietary-related viruses in the gut virome of marmosets hybrids (Callithrix jacchus x Callithrix penicillata), Brazil.

Viral metagenomics has contributed enormously to the characterization of a wide range of viruses infecting animals of all phyla in the last decades. Among Neotropical primates, especially those introduced, knowledge about viral diversity remains poorly studied. Therefore, using metagenomics based on virus enrichment, we explored the viral microbiota present in the feces of introduced common marmosets (Callithrix sp.) in three locations from the Silva Jardim region in the State of Rio de Janeiro, Brazil. Fecal samples were collected from nine marmosets, pooled into three sample pools, and sequenced on Illumina MiSeq platform. Sequence reads were analyzed using a viral metagenomic analysis pipeline and two novel insect viruses belonging to the Parvoviridae and Baculoviridae families were identified. The complete genome of a densovirus (Parvoviridae family) of 5,309 nucleotides (nt) was obtained. The NS1 and VP1 proteins share lower than 32% sequence identity with the corresponding proteins of known members of the subfamily Densovirinae. Phylogenetic analysis suggests that this virus represents a new genus, provisionally named Afoambidensovirus due to its discovery in the Brazilian Atlantic Forest. The novel species received the name Afoambidensovirus incertum 1. The complete circular genome of a baculovirus of 107,191 nt was also obtained, showing 60.8% sequence identity with the most closely related member of the Baculoviridae family. Phylogenetic analysis suggests that this virus represents a new species in the Betabaculovirus genus, provisionally named Betabaculovirus incertum 1. In addition, sequences from several families of arthropods in the three pools evaluated were characterized (contigs ranging from 244 to 6,750 nt), corroborating the presence of possible insect hosts with which these new viruses may be associated. Our study expands the knowledge about two viral families known to infect insects, an important component of the marmosets' diet. This identification in hosts' feces samples demonstrates one of the many uses of this type of data and could serve as a basis for future research characterizing viruses in wildlife using noninvasive samples.

Multimodal analyses of a non-human primate model harboring mutant amyloid precursor protein transgenes driven by the human EF1α promoter.

Alzheimer's disease (AD) is the leading cause of dementia which afflicts tens of millions of people worldwide. Despite many scientific progresses to dissect the AD's molecular basis from studies on various mouse models, it has been suffered from evolutionary species differences. Here, we report generation of a non-human primate (NHP), common marmoset model ubiquitously expressing Amyloid-beta precursor protein (APP) transgenes with the Swedish (KM670/671NL) and Indiana (V717F) mutations. The transgene integration of generated two transgenic marmosets (TG1&TG2) was thoroughly investigated by genomic PCR, whole-genome sequencing, and fluorescence in situ hybridization. By reprogramming, we confirmed the validity of transgene expression in induced neurons in vitro. Moreover, we discovered structural changes in specific brain regions of transgenic marmosets by magnetic resonance imaging analysis, including in the entorhinal cortex and hippocampus. In immunohistochemistry, we detected increased Aß plaque-like structures in TG1 brain at 7 years old, although evident neuronal loss or glial inflammation was not observed. Thus, this study summarizes our attempt to establish an NHP AD model. Although the transgenesis approach alone seemed not sufficient to fully recapitulate AD in NHPs, it may be beneficial for drug development and further disease modeling by combination with other genetically engineered models and disease-inducing approaches.

Integrated analysis of long-noncoding RNA and circular RNA expression in Peste-des-Petits-Ruminants Virus (PPRV) infected marmoset B lymphocyte (B95a) cells.

Peste-des-Petits-Ruminants (PPR) or goat plague is an important viral disease of sheep and goats caused by the small ruminant morbilli virus or PPR virus (PPRV). Long non coding RNAs (lncRNA) and circular RNAs (circRNA) play a pivotal role in several biological processes including regulation of virus-host interactions. The present study explored the expression of lncRNA, circRNA and their functions in PPRV infected B-lymphocyte (B95a) cells. The results revealed a total of 4531 lncRNA and 2348 circRNA expression in both mock and PPRV infected samples. Analysis of differentially expressed (DE) RNA identified 123 DE-lncRNA and 39 DE-circRNA as significantly dysregulated. Functional analysis of cis-target genes of DE-lncRNA indicated activation of TCF dependent WNT signaling and PKN1 stimulated transcription process. Interactions (sponging) of microRNA (miRNA) revealed 344 DE-lncRNA-miRNA and 93 DE-circRNA-miRNA pairs. The competing endogenous RNA (ceRNA) network of lncRNA/circRNA-miRNA-mRNA in PPRV infected B95a cells was represented by 69 ceRNA pairs. We validated the DE-circRNA by targeted amplification and sequencing of back spliced junctions (BSJs). The present study revealed a profile of lncRNA, circRNA and their potential ceRNA network in PPRV infection. The results provide insight for better understanding of PPRV-host interactions.

Intraintestinal Analysis of the Functional Activity of Microbiomes and Its Application to the Common Marmoset Intestine.

The intestinal microbiome is closely related to host health, and metatranscriptomic analysis can be used to assess the functional activity of microbiomes by quantifying microbial gene expression levels, helping elucidate the interactions between the microbiome and the environment. However, the functional changes in the microbiome along the host intestinal tract remain unknown, and previous analytical methods have limitations, such as potentially overlooking unknown genes due to dependence on existing databases. The objective of this study is to develop a computational pipeline combined with next-generation sequencing for spatial covariation analysis of the functional activity of microbiomes at multiple intestinal sites (biogeographic locations) within the same individual. This method reconstructs a reference metagenomic sequence across multiple intestinal sites and integrates the metagenome and metatranscriptome, allowing the gene expression levels of the microbiome, including unknown bacterial genes, to be compared among multiple sites. When this method was applied to metatranscriptomic analysis in the intestinal tract of common marmosets, a New World monkey, the reconstructed metagenome covered most of the expressed genes and revealed that the differences in microbial gene expression among the cecum, transverse colon, and feces were more dynamic and sensitive to environmental shifts than the abundances of the genes. In addition, metatranscriptomic profiling at three intestinal sites of the same individual enabled covariation analysis incorporating spatial relevance, accurately predicting the function of a total of 10,856 unknown genes. Our findings demonstrate that our proposed analytical method captures functional changes in microbiomes at the gene resolution level. IMPORTANCE We developed an analysis method that integrates metagenomes and metatranscriptomes from multiple intestinal sites to elucidate how microbial function varies along the intestinal tract. This method enables spatial covariation analysis of the functional activity of microbiomes and accurate identification of gene expression changes among intestinal sites, including changes in the expression of unknown bacterial genes. Moreover, we applied this method to the investigation of the common marmoset intestine, which is anatomically and pharmacologically similar to that of humans. Our findings indicate the expression pattern of the microbiome varies in response to changes in the internal environment along the intestinal tract, and this microbial change may affect the intestinal environment.

Region-specific microRNA alterations in marmosets carrying SLC6A4 polymorphisms are associated with anxiety-like behavior.

BACKGROUND: Psychiatric diseases such as depression and anxiety are multifactorial conditions, highly prevalent in western societies. Human studies have identified a number of high-risk genetic variants for these diseases. Among them, polymorphisms in the promoter region of the serotonin transporter gene (SLC6A4) have attracted much attention. However, due to the paucity of experimental models, molecular alterations induced by these genetic variants and how they correlate to behavioral deficits have not been examined. In this regard, marmosets have emerged as a powerful model in translational neuroscience to investigate molecular underpinnings of complex behaviors. METHODS: Here, we took advantage of naturally occurring genetic polymorphisms in marmoset SLC6A4 gene that have been linked to anxiety-like behaviors. Using FACS-sorting, we profiled microRNA contents in different brain regions of genotyped and behaviorally-phenotyped marmosets. FINDINGS: We revealed that marmosets bearing different SLC6A4 variants exhibit distinct microRNAs signatures in a region of the prefrontal cortex whose activity has been consistently altered in patients with depression/anxiety. We also identified Deleted in Colorectal Cancer (DCC), a gene previously linked to these diseases, as a downstream target of the differently expressed microRNAs. Significantly, we showed that levels of both microRNAs and DCC in this region were highly correlated to anxiety-like behaviors. INTERPRETATION: Our findings establish links between genetic variants, molecular modifications in specific cortical regions and complex behavioral responses, providing new insights into gene-behavior relationships underlying human psychopathology. FUNDING: This work was supported by France National Agency, NRJ Foundation, Celphedia and Fondation de France as well as the Wellcome Trust.

Early development of the cochlea of the common marmoset, a non-human primate model.

BACKGROUND: Fine-tuned cochlear development is essential for hearing. Owing to the difficulty in using early human fetal samples, most of our knowledge regarding cochlear development has been obtained from rodents. However, several inter-species differences in cochlear development between rodents and humans have been reported. To bridge these differences, we investigated early otic development of a non-human primate model animal, the common marmoset (Callithrix jacchus). METHODS: We examined 20 genes involved in early cochlear development and described the critical developmental steps for morphogenesis, which have been reported to vary between rodents and marmosets. RESULTS: The results revealed that several critical genes involved in prosensory epithelium specifications showed higher inter-species differences, suggesting that the molecular process for hair cell lineage acquisition in primates differs considerably from that of rodents. We also observed that the tempo of cochlear development was three times slower in the primate than in rodents. CONCLUSIONS: Our data provide new insights into early cochlear development in primates and humans and imply that the procedures used for manipulating rodent cochlear sensory cells cannot be directly used for the research of primate cells due to the intrinsic inter-species differences in the cell fate determination program.

The gut microbiome of exudivorous marmosets in the wild and captivity.

Mammalian captive dietary specialists like folivores are prone to gastrointestinal distress and primate dietary specialists suffer the greatest gut microbiome diversity losses in captivity compared to the wild. Marmosets represent another group of dietary specialists, exudivores that eat plant exudates, but whose microbiome remains relatively less studied. The common occurrence of gastrointestinal distress in captive marmosets prompted us to study the Callithrix gut microbiome composition and predictive function through bacterial 16S ribosomal RNA V4 region sequencing. We sampled 59 wild and captive Callithrix across four species and their hybrids. Host environment had a stronger effect on the gut microbiome than host taxon. Wild Callithrix gut microbiomes were enriched for Bifidobacterium, which process host-indigestible carbohydrates. Captive marmoset guts were enriched for Enterobacteriaceae, a family containing pathogenic bacteria. While gut microbiome function was similar across marmosets, Enterobacteriaceae seem to carry out most functional activities in captive host guts. More diverse bacterial taxa seem to perform gut functions in wild marmosets, with Bifidobacterium being important for carbohydrate metabolism. Captive marmosets showed gut microbiome composition aspects seen in human gastrointestinal diseases. Thus, captivity may perturb the exudivore gut microbiome, which raises implications for captive exudivore welfare and calls for husbandry modifications.

Uniparental genetic markers to investigate hybridization in wild-born marmosets with a mixed phenotype among Callithrix aurita and invasive species.

The native marmoset of the Southeastern Atlantic Forest in Brazil is among the 25 most endangered primates of the world. Hybridization with alien species is one of its main threats registered since the early 2000s based on phenotype, so far, without genetic confirmation. Using uniparental molecular markers, we analyzed 18 putative hybrids, captured from 2004 to 2013 in different localities of the Atlantic Forest. A nine base pair deletion in the SRY gene of C. aurita was used to investigate paternal ancestry. Maternal ancestry was assessed by DNA sequencing of ca. 455 bp from the COX2 gene. Hybridization was confirmed for 16 out of the 18 marmosets since they inherited COX2 haplotypes of the alien C. penicillata or C. jacchus and the SRY deletion specific to C. aurita. Two individuals inherited both parental lineages of C. aurita, which is probably related to backcrossing or hybrid interbreeding. The direction of hybridization of females with the matrilineal lineage of invasive species with males descending from the native lineage was predominant in our sampling. This is the first time that hybridization between C. aurita and invasive species has been confirmed through genetic analysis.

AAV capsid variants with brain-wide transgene expression and decreased liver targeting after intravenous delivery in mouse and marmoset.

Genetic intervention is increasingly being explored as a therapeutic option for debilitating disorders of the central nervous system. The safety and efficacy of gene therapies rely upon expressing a transgene in affected cells while minimizing off-target expression. Here we show organ-specific targeting of adeno-associated virus (AAV) capsids after intravenous delivery, which we achieved by employing a Cre-transgenic-based screening platform and sequential engineering of AAV-PHP.eB between the surface-exposed AA452 and AA460 of VP3. From this selection, we identified capsid variants that were enriched in the brain and targeted away from the liver in C57BL/6J mice. This tropism extends to marmoset (Callithrix jacchus), enabling robust, non-invasive gene delivery to the marmoset brain after intravenous administration. Notably, the capsids identified result in distinct transgene expression profiles within the brain, with one exhibiting high specificity to neurons. The ability to cross the blood-brain barrier with neuronal specificity in rodents and non-human primates enables new avenues for basic research and therapeutic possibilities unattainable with naturally occurring serotypes.

Efficient Induction of Primate iPS Cells Using a Combination of RNA Transfection and Chemical Compounds.

Non-human primate induced pluripotent cells (iPS cells) are useful for preclinical studies of iPS cell-based therapies and the research of primate developments. Since the initial report of iPS cells in 2006, various iPS cell induction methods have been reported. Here, we describe an efficient method for inducing iPS cells using a combination of RNA transfection and chemical compounds without using transgenes. Many kinds of marmoset cells, including difficult-to-reprogram cells, can be converted into iPS cells using this combinatorial method. Furthermore, this method can be applied to other primates, including humans.